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Abstract 

By applying Darboux-Crum transformations to a Lax pair formulation of the Korteweg-de 
Vries (KdV) equation, we construct new sets of multi-soliton solutions to it as well as to the 
modified Korteweg-de Vries (mKdV) equation. The obtained solutions exhibit a chiral asym¬ 
metry in propagation of different types of defects in crystalline backgrounds. We show that 
the KdV solitons of pulse and compression modulation types, which support bound states in, 
respectively, semi-infinite and finite forbidden bands in the spectrum of the perturbed quantum 
one-gap Lame system, propagate in opposite directions with respect to the asymptotically peri¬ 
odic background. A similar but more complicated picture also appears for multi-kink-antikink 
mKdV solitons that propagate with a privileged direction over the topologically trivial or topo¬ 
logically nontrivial crystalline background in dependence on position of energy levels of trapped 
bound states in spectral gaps of the associated Dirac system. Exotic N = 4 nonlinear super- 
symmetric structure incorporating Lax-Novikov integrals of a pair of perturbed Lame systems 
is shown to underlie the Miura-Darboux-Crum construction. It unifies the KdV and mKdV 
solutions, detects the defects and distinguishes their types, and identifies the types of crystalline 
backgrounds. 


1 Introduction 

Nonlinear integrable systems play an important role in a variety of areas of physics and its applica¬ 
tions [31313II El 13 HE]- a Lax pair formulation in particular cases of the Korteweg-de Vries (KdV) 
and modified Korteweg-de Vries (mKdV) equations relates their soliton, kink, kink-antikink, and 
crystalline type solutions with reflectionless and finite-gap systems of quantum mechanics mmm- 
This formulation for the KdV and mKdV systems includes, respectively, the stationary Schrodinger 
and Dirac equations (the latter with a scalar potential) as the eigenstate equation. 

The inverse scattering method allows to construct soliton solutions for the KdV equation in 
analytical form for the asymptotically free, homogeneous boundary conditions. On the other hand, 
such analytic solutions can also be constructed by exploiting the covariance of the Lax pair under 
Darboux transformations and their generalization in the form of Darboux-Crum transformations. 
From the spectral point of view, these transformations correspond to adding soliton defects by 
means of introducing bound states into the spectrum of the initial potential [ID1 . 

Recently in |12| . by employing Darboux-Crum transformations, we constructed a generalization 
of reflectionless potentials by introducing soliton defects into the periodic (crystalline) background 
of the one-gap Lame system 

Ui. nrn e(x) = 2k 2 sv? (x\k) + const. (1-1) 
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As a result, we obtained one-gap potentials for the Schrodinger system with an arbitrary number 
of bound states in the lower and intermediate forbidden bands, which are trapped by the soliton 
defects. 

Using a relation between Darboux transformations and supersymmetric quantum mechanics, it 
was also possible to construct finite-gap Dirac (Bogoliubov-de Gennes) Hamiltonian operators, in 
which scalar potentials carry perturbations of the kink and kink-antikink type introduced into the 
crystalline backgrounds of different nature. 

One such periodic background corresponds to the superpotential of the form da m m, \w\ 


V KC (x) 


— log dn (x \k) = — k 2 
ax 


cn (x | k) sn {x \k) 
dn (x \ k) 


( 1 . 2 ) 


This finite-gap Dirac scalar potential appeared as a solution in the Gross-Neveu model with a dis¬ 
crete chiral symmetry mm, and was identified there as the kink crystal. In the QCD framework, 
the kink crystal can be considered as the phase corresponding to the crystalline color superconduc¬ 
tor m, that is the QCD analogue of the Larkin-Ovchinnikov-Fulde-Ferrell phase in the context 
of electron superconductivity mm- The kink crystal type solution (11.21) also found some appli¬ 
cations in the physics of conducting polymers [221 [23 ]. In the construction from m, the Dirac 
Hamiltonian operator with the perturbed scalar potential of the form m possesses a spectrum 
having a central allowed band between the conduction band and the Dirac sea, and the bound 
states appear there symmetrically in both forbidden bands (gaps). 

There is yet another possibility for asymptotic behaviour in the form of the kink-antikink, V KA , 
(or, antikink-kink, V AK ) crystal PI, 

V KA (x) = Z(x + /3\k)-Z(x\k)-z((3\k), V AK (x) = -V KA (x ), (1.3) 


which was found as a solution in the Gross-Neveu model with a bare mass term [24] . It allowed 
to obtain scalar potentials with two permitted bands between the Dirac sea and conduction band, 
with arbitrary number of bound states in central and/or other two forbidden bands. Here z(/3\k) 
is a constant given in terms of the Jacobi Zeta and elliptic functions, z(/3\k) = Z(j3 + iK 7 \k) + 
i = Z{P\k) + cn(P\k)dn{P\k)/sn (0|fc). An additional change of the global topology of the 
Dirac scalar potential provided us with reflectionless kink-antikink perturbations over the kink-type 
crystalline background. 

A natural question we address in this paper is the following: how the dependence on time can be 
‘reconstructed’ for the finite-gap Schrodinger and Dirac crystalline potentials with defects, found 
in m, so that they will be solutions for the KdV and mKdV equations? We refer here to the 
evolution parameter in the corresponding nonlinear integrable systems; from the perspective of the 
eigenvalue problems, it is associated with the isospectral deformation of the potentials. Answering 
this question, we reveal a phenomenon of a chiral asymmetry in propagation of different types of 
soliton defects in crystalline backgrounds, which, as we believe, could be interesting, particularly 
from the point of view of applications. 

The paper is organized as follows. In the next Sectional that is of a very brief review character, 
we consider the Lax pair formulation for the KdV equation, discuss its Darboux covariance, and 
apply the Darboux and Darboux-Crum transformations to obtain multi-soliton solutions for the 
KdV equation. Section [3] is devoted to the construction of solutions to the auxiliary problem 
corresponding to the Lax pair with the one-gap Lame potential taken as the stationary, periodic 
cnoidal solution to the KdV equation. On the basis of the solutions to the auxiliary problem, in 
Section[]]we construct solutions to the KdV equation by employing Darboux-Crum transformations. 
In this way we obtain the KdV solutions that describe propagation of solitons of the potential well 
(pulse) and compression modulation types in the background of asymptotically periodic cnoidal 
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wave. We also discuss there the issue of velocities of the defects. In Section [5l we first show 
how the Miura and Darboux-Crum transformations together with the Galilean symmetry of the 
KdV equation can be employed for the construction of solutions to the mKdV equation on the 
basis of the KdV solutions. Then we discuss the exotic N = 4 nonlinear supersymmetric structure 
incorporating Lax-Novikov integrals of a pair of perturbed Lame systems. This supersymmetry 
underlies the Miura-Darboux-Crum construction and unifies the KdV and mKdV solutions. In 
Section[6]we first discuss briefly the asymptotically free mKdV solutions corresponding to the multi- 
kink-antikink solitons propagating over the kink or kink-antikink background. Then we construct a 
much more rich set of multi-kink-antikink soliton solutions for the mKdV equation which propagate 
over the topologically trivial or topologically nontrivial crystalline backgrounds. The last Section 
[7] is devoted to the discussion and outlook. 

2 Lax pair for the KdV equation, Darboux transformations, and 
multi-soliton solutions 

2.1 Auxiliary spectral problem for the KdV equation 

Consider a linear system 

L<j> = A0, ^ = P<j) (2.1) 

for a function 0 = (f>(x,t, A). It is assumed that L and P are some (in general, matrix) differential 
operators in space coordinate i£K with coefficients that can also depend on evolution parameter 
t. If the evolution in t generated by P is isospectral, dX/dt = 0, then the condition of consistency 
for the system (12.ip reduces to the Lax equation 

r)T 

^ = [P,L]. (2.2) 

For the choice of the Lax pair in the form of differential operators 

L = ~d 2 x + u, (2.3) 

P = —45)! + 6ud x + 3 u x , (2.4) 

equation (12.2j) reduces to the Korteweg-de Vries equation for the scalar field u = u(x,t), 


U £ — GuU x ^xxx ■ 


2.2 Darboux covariance of the KdV equation 

The system of equations corresponding to the Lax pair (12.31) . (12.41) . 

(—d 2 + u)'ff(x,t,X) = AT(x, t, A) , 

5 q 

—T(x,t, A) = (-45)! + 6ud x + 3u x )^(x, t, A), 
is covariant under Darboux transformations Coj 


u(x,t ) 
V(x,t, A) 


ui{x,t) = 
^i(x,t, A) 


u(x,t) - 2(log'F(x,f, Ai)) xx , 
_ W($(x,t, Xi),^{x,t, A)) 
V(x,t, Ai) 


(2.5) 


( 2 . 6 ) 

(2.7) 


( 2 . 8 ) 

(2.9) 
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where W is the Wronskian, W (/, g ) = fg x — f x g■ This follows from the observation that if Eqs. (12.61) 
and (12.71) are fulfilled for T(x, t, A), and 11/(x, t, Ai) satisfies the same equations with A changed for Ai, 
then Ti(x, t, A) obeys the equations (12.61) . (12.71) with u(x, t) changed for u\(x, t). As a consequence, 
if u(x,t) is a solution of the KdV equation (|2.5|) . then u\(x,t) obeys the same equation. 

This result can be extended for a finite sequence of consecutive Darboux transformations, 


u(x,t ) 
V(x,t, A) 


Um(x,t ) = u(x,t) - 2(log W (^(x, t, Xi), . . . ,T(x,f, \m)))xx , 
TT(T(x,t,Ai),...,^(x,t,A m ),T(x,t,A)) 
m[X ' ’ j lT(^(x,t,A 1 ),...,T(x,t,A m )) 


( 2 . 10 ) 

( 2 . 11 ) 


that is known as the Darboux-Crum transformation of order m |10j . For the sake of simplicity of 
notations, the dependence of u m (x, t ) as well as of T m (x, t, A) on the spectral parameters Ai,..., A m 
is not shown explicitly. 

One can define a first order differential operator A\ in terms of the function \k(x,t, Ai) that 
underlies the Darboux transformation construction (12.81) . (12.91) . 


A\ = T(x, t, \i)d x 1 =d x ~Vi (x,f), Fi(x,t) = (logT(x,t, Ai)) x . (2.12) 

v(x,t, Ai) 

The operator A\ and its Hermitian conjugate, A \, intertwine the Schrodinger operators L = —d% + 
u(x, t ) and L\ = —d\ + u\(x, t ), 


A\L — LiAi, 


a\l x = la\, 


and factorize them, 

= L — Ai, = L\ — Ai. 

The two relations in ((2.14)1 are equivalent to Riccati equations for the function Vj, 


(2.13) 

(2.14) 


V? + V lx = u- X u V? - V lx = Ul - Ai. 


(2.15) 


For the successive Crum-Darboux transformations of orders m — 1 and m, the corresponding 
Schrodinger operators L m _i and L m are related by means of the first order differential operator 

-A m — (A m _xT(x, t, A m))d x tt ~r~, — r 77 d x V? n (x, t), m 1,... , (2.16) 

(A m -lW(x,t, Am)) 

and its conjugate, where A n = A n A n _i ... Ai, Ao = 1. Then, as a generalization of (12.13)1 and 
(|^A4p. we have A rn L rri _\ — Z/ m A m , A m Tm — T m—iA-m , and A m A m — L m _\ A A m A m — 

Lm - Am, where L 0 = L. 

Superpotential V m can be presented in terms of the Wronskians p5] . 

V m (x,t) = n m _i(x,t) - D m (x,t), = - (log W (4'(x, t, Ai),... ,^(x,t, A m ))) x , (2.17) 

where we assume Do = 0 and W(^(x, t, Ai)) = T(x,t, Ai). It satisfies the Riccati equations 

Vm(x,t ) + (V m (x,t)) x = u m -i(x,t) - A m, V^(x,t) - (V m (x,t)) x = u m (x,t ) - A m . (2.18) 
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2.3 Multi-soliton solutions of the KdV equation 


The described picture with Darboux-Crum transformations can be illustrated by the construction 
of multi-soliton solutions for the KdV equation. 

The trivial solution of the KdV equation is uo = 0. In this case Lq = — corresponds to 
the Schrodinger operator for a free particle system, and the evolution operator (El reduces to 
Pq = — 4^3-. The system (12.11) takes then the form 


<9 2 T _ <9T d 3 T 

-= AT -= —4- 

dx 2 1 dt dx 3 


(2.19) 


Acting on both sides of the first equation by 4d x and summing up the result with the second 
equation, we obtain = 4Ad x 4/. Therefore, T(x, t, A) = T(x + 4A t, A). For A = 0, A = — k 2 < 0, 
and A = k 2 > 0, the pairs of linearly independent solutions of the system (12.191) can be chosen in 
the form 

= 0) = {1, x}, (2.20) 

T(x,f, A = — k 2 ) = {cosh X~, sinh A - }, T(x, t, A = k 2 ) = {cos X + , sinA + }, (2.21) 

where = k(x — xq =f 4 n 2 t). To apply the Darboux-Crum transformation to a trivial solution 
uo = 0, we have different possibilities for choosing wave functions from the sets ()2.20l) . (12.211) . The 
choice 4>(x,t, 0) = x gives rise by means of the first order Darboux transformation (12.81) to the 
simplest but singular, time-independent solution of the KdV equation, u\{x) = 2/a; 2 . The non¬ 
singular solutions for the KdV equation are generated by choosing appropriately the eigenstates 
with eigenvalues A < 0, 

d 2 

u n (x,t) = “2-^2 log IV(coshX 1 ,sinhX 2 ..., f{X n )), X- = Kj (x - x 0j - 4n 2 t), (2.22) 


where the last argument in the Wronskian is f(X~) = sinh X~ if n is even, n = 21, and f(X~) = 
cosh X~ for odd n = 21 + 1; xqj are the translation (phase) parameters, and the scale parameters 
Kj have to obey the inequalities 0 < ki < K2 < k^ < ... < K n . Functions (12.221) correspond to 
n-soliton solutions of the KdV equation. The case n = 2 is illustrated in Fig. [T] When solitons 
in the solution (12.221) are well separated, the propagation to the right of the j-th soliton can be 
characterized by the speed Vj = 4k 2 and amplitude 2k 2 . 


3 Spectral problem with the cnoidal background 

The simplest stationary periodic solution to the KdV equation (12.51) can be presented in a form 

2 

u{x) = UQfl{x) = 2&- 2 /r 2 sn 2 (jix \k) — - (l + k 2 ) /a 2 , (3.1) 

o 

where sn (u\k) is the Jacobi elliptic function, whose real and imaginary periods depend on the 
modular parameter 0 < A; < 1, //>0isa free (scale) parameter, and the sense of the indices in uo,o 
will be clarified below. Due to the f-independence of solution (13.11) . Lax equation (12.21) reduces to 
the condition of commutativity of the corresponding operators (12.31) and (12.41) constructed on the 
basis of (13.11) . 

[L,P]= 0. (3.2) 

Relation (13.21) guarantees the existence of the common basis for the operators L and P. We then 
look for the solutions of the system of equations ()2.6I) . (12.71) in the form 

T(x, t, A) = T(x, a) exp(n(a)t ), (3-3) 
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Figure 1: The KdV two-soliton solution with K\ = 1, k 2 = 1.4, and xq i = 0, i = 1,2, is shown by 
a continuous line. Propagation of a one-soliton solution with k,\ = 1 and xoi = 0 is depicted by a 
dashed line. Initial phases are chosen so that the two-soliton and one soliton solutions at t = 0 are 
symmetric with respect to the point x = 0. On the right, the spectrum of the associated two-soliton 
Schrodinger operator L 2 is shown. The continuous green line corresponds to a doubly degenerate 
semi-infinite continuos part of the spectrum with eigenstates ^^(x, t) = A 2 e zizlX+ ^ x ’ t ' ,K ’ Xo \ while 
the filled circle indicates a non-degenerate edge-state described by the eigenfunction ipo(x, t ) = A 2 I. 
Dashed red line corresponds to non-physical semi-infinite part of the spectrum, inside which blue 
squares indicate energies of the two bound states trapped by solitons and described by eigenfunc¬ 
tions i/>i(x,t ) = A 2 sinhXf and ip 2 (x,t) = A 2 CoshX^~. 


where <h(x, a) is a common eigenstate of L and P , 3?(x, a) = A(ct)4>(a;, a), P<h(x, a) = 7 r(o!)$(x, a). 

The sought for state <h(x, a) is [12] [16] [26j 


= H (fix + a\k) ux z(g\k) 
0 (fix\k) 


(3.4) 


where H, 0 and Z are Jacobi’s Eta, Theta and Zeta functions, while the corresponding eigenvalues 
are 

X(a\k) = fi 2 (dn 2 (a|fc) — ^(1 + k' 2 )) , (3.5) 

n(a\k) = — 4/c 2 /r 3 sn (a \k) cn (a |fc) dn (a \k) . (3-6) 

Notice that A(— a) = A(a), 7 r(—a) = — 7r(a), 7r(a) = 2; an< 4 

n\a) = ~16(A(a) -E 0 )(A(a) -Ei)(A(a) -E 2 ), (3.7) 

where 

E 0 = -^(l + fc'V, = |(1 - 2A 2 )/X 2 , E 2 = l -(l + k 2 )Li 2 , (3.8) 

and k' = \/l — k 2 is the complementary modular parameter. Relations (13.21) and (13.71) correspond 
to the Burchnall-Chaundy theorem m , and reflect the one-gap nature of the potential 0 mm, 
see below. 
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Taking into account the double periodicity of the Jacobi sn, cn and dn functions, without loss 
of generality one can suppose that a dimensionless parameter a takes values in the rectangle in 
the complex plane with vertices irfl 0, K, K + zK 7 and zK 7 . On the border of this fundamental 
a-rectangle, eigenvalues A(ct) of the Schrodinger operator L take real values, while the eigenvalue 
ir(a) of the Lax operator P takes nonzero real values on the horizontal borders of this rectangle 
where a = (3 + , 0 < /3 + < K, and a = f3~ + zK 7 , 0 < < K. On the other hand, on both vertical 

borders of the rectangle, a = ry" 1- , 0 < 7 + < K 7 , and a = K + ry _ , 0 < 7 ^ < K 7 , 7 r(a) takes 
pure imaginary values, turning into zero at the vertices 0, K and K + zK 7 . The vertex zK 7 is the 
third order pole of 7 r(a) and the second order pole of the function A(a). Moreover, corresponding 
eigenfunctions (13.41) are unbounded functions on the horizontal edges, but they are bounded on the 
vertical borders of the rectangle. Summarizing, in correspondence with the specified properties, 
the solutions (|3.3I) of the auxiliary spectral problem for the KdV equation are bounded on the 
vertical borders of the rectangle, except the vertex point a = zK 7 , and they are unbounded on the 
horizontal borders except the vertices a = 0, K, K + zK 7 . 

The Schrodinger operator L with potential (|3.ip describes quantum mechanical periodic one- 
gap Lame system, for which the upper horizontal line a = f5~ +zK' corresponds to the semi-infinite 
forbidden band with energy values —00 < E(a ) = A (a) < Eq, and the lower horizontal line a = /3 + 
corresponds to a spectral gap with E\ < E(a) < Eo. The corresponding bounds Eq, E\ and E 2 , 
defined in (I3.8[) . satisfy the relations Eq < 0, £2 > 0, while E\ can take positive or negative values 
in dependence on the value of the modular parameter k. The vertical borders of the rectangular 
correspond to the allowed valence, a = K + zy - , E$ < E(a) < E \, and conduction, a = Z 7 + , 
E 2 < E(a) < 00 , bands. With taking into account relations A(— a) = A(a), 7r(— a) = —zr(a), the 
common eigenstates of L and P with opposite values for 7 r(a) are obtained from the described 
solutions with a restricted to the borders of the indicated rectangle by changing a —>■ —a. 

Let us stress that the solutions (13.31) to the system of equations (12.61) and (12.71) have a factorizable 
dependence on the evolution variable t due to the f-independence of the corresponding KdV solution 
(13.11) and the associated with it Lax pair (j2.3|) . (|2.4D . 

4 Multi-soliton defects in the crystalline background 

One can construct three types of Darboux-Crum transformations based on the solutions (13. 3 p to 
the auxiliary problem for the KdV equation. They will provide us with new soliton solutions of 
the KdV equation that propagate over the stationary crystalline background (13.11) by deforming it. 
These are: 

i) soliton defects of the potential well (pulse) type; 

ii) soliton defects of the compression modulations nature; 

iii) the mixed case in which both types of the solitons are present. 

Soliton defects of the type i) are generated by Darboux-Crum transformations based on unbounded 
(non-physical) states from the lower semi-infinite forbidden band in the spectrum of the one-gap 
Lame system. The construction of soliton defects of the type ii) requires the use of states from the 
gap of the quantum Lame system. The employment of both types of non-physical states from the 
spectrum of the Lame system generates solutions corresponding to the mixed case iii). Below we 
describe these three types of solutions. 

1 Here, K = K(fc) is a complete elliptic integral of the first kind and K' = K(fc'). 
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4.1 Potential well (pulse) type solitons 

If we choose common eigenstates (13.31) of L and P in the form of real-valued functions, then the new 
solutions generated by Darboux-Crum transformations (12.101) also will be real. The eigenfunctions 
<h(x,a = /3~ + iK'), f3~ € (0, K), of L in (13.4p correspond to non-physical (unbounded) states with 
negative energies from the lower forbidden band of one-gap Lame system. They are real-valued 
functions modulo a global phase factor mm- Omitting a phase factor, we obtain a common 
real-valued eigenstates of L and P, 


F(x,t,/3 ) 


9 ^© X ( + Lx\k) k ^ 6XP l fc ) + 7r (/ 3 +iK'|fc)i), 


(4.1) 


where 


z (r \k) = z (r + iK'\k)+i^- = z (r i*o + 


cn(/3 \k)dn(/3 \k) 


2K 


sn (/3 | k ) 


> 0 , 


7t(/3 + iK.'\k) = 4/r' 


3 cn (ft |fc)dn(/3 |fc) 
sn 3 (ft - |fc) 


> 0 . 


(4.2) 

(4.3) 


The superpotential (11.31) is given in terms of (|4.1I) by V KA (x) = (log F(x, t, /3)) x with ft = ft - and 
H = 1, that reduces to (11.21) in the limit case ft - = K. On the basis of the functions (14.11) . we 
construct Darboux-Crum transformation of the form 


u 0 ,i ( x >*) = “ 0 , 0 ( 2 ;) - 2(logIT(J r + (l), J 7 _(2), ...,F(_iy+i(l))) xx , (4.4) 

where 

F±{j) = CjF(x,t,Pj)±-^=F(x,t,-Pj), j = (4.5) 

are combinations of the two linear independent solutions (14.11) . Minus index in parameter coefficients 
CJ indicates that (14.11) are real eigenstates from the lower forbidden band of the one-gap Lame 
system. The choice 

K > Pi > P 2 > ••• > /3 Z _ > 0 , 0 < CJ < 00 , (4.6) 

guarantees the non-singular nature of the solution uoj(x, t ) for the KdV equation [12]. This solution 
describes l solitons of the potential well type that propagate to the right in the stationary periodic 
background by deforming it. For large negative and positive values of t, pulses are well separated, 
and each corresponding potential well supports a bound state with negative energy given by (13.51) 
with a = Oij = /3j + iK 7 . A deeper potential well soliton defect in asymptotically periodic back¬ 
ground propagates faster, and supports the bound state with lower energy, see Section 14.41 below. 
Figures [2] and [3] show such travelling soliton defects for the simplest cases of l = 1 and 1 = 2. 

In described solutions, an arbitrary number of solitons can be eliminated by taking limits of the 
type CJ —>• 0 or CJ —>• 00 that corresponds to sending the j-th soliton to minus or plus infinity. 
This provokes a global phase shift x —>• x — -Jr for C~ —>• 0, or x ^ x + -J. w hen CJ —>• 00 , in the 

J 3 &3 

crystalline background as well as in the remaining solitons, and additional change in the parameters 
CJ,, j' / j, which depend on /3j, and /3j, see refs. mm- 

The indicated phase shifts of the crystalline background can be understood in the following 
simple way. Consider the first order Darboux transformation based on the function F{x,fi~) 
from Eq. (14.11) and applied to the stationary periodic KdV solution (13.11) . Since log F(x,/3~) = 
log 0 (fix + /3~ | k ) — log 0 (fix | k ) — (/rxz(/3“ \k) — 7r(/3“ + iK'\k)t ), and — 2 (log 0 (fix\k)) xx = 
“ 0 , 0 ( 2 ;), we see that the Darboux transformation generated by F(x,/3~) transforms 110 , 0 ( 2 ;) into 
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U 0 ,l(x,t) 


a (^0,l) 



Figure 2: The KdV solution with a pulse soliton propagating over the stationary crystalline back¬ 
ground ()3. 1 1) shown by dashed line. Here /x = 1, k = 0.6, = 1.2, Cj - = 1. Vertical line on the 

right illustrates the spectrum of the perturbed Lame system with lower and upper forbidden bands 
shown in dashed red. Filled semicircles indicate non-degenerate energy values at the edges of the 
valence and conduction bands. In the semi-infinite lower forbidden band, there is a bound state 
trapped by the soliton defect, that is shown by a blue square. Parameters in (j4.5)l are chosen so 
that the solution at t = 0 is symmetric under the space reflection x —>• —x. 


UQfl{x + /3 /fx), while the transformation based on F(x,—/3 ) produces the displacement x — » 
x — /H- When in the Darboux-Crum transformation (14.41) we take the limit Cj —» 0, in the 

function (J4.5I) the term with F(x, — /3~) survives. Due to relation z(—/3~\k) < 0, this function 
exponentially increases in the region x —> +oo, and this provokes the displacement x —>• x — /3 _ //x 
of the asymptotically periodic background of the solution in that region. Analogously, for the 
limit Cj —>• oo, the first term with F(x,/3j) survives in the function (14.51) . and the asymptotically 
periodic background in the solution will be displaced by the shift x —> x + /3 _ //x in the region to 
the left {x —>• —oo) from the rest of the surviving solitons. 

4.2 Solitons of the compression modulation type 

The states (13.31) with a = ±/3 + , 0 < /3 + < K, which we denote here as 4>(a:, t, ±/3 + ), 

4>(x, t, /3 + ) = H ex P ( “ W Z(/3 + | k) + ir{/3 + \k)t s ) , (4.7) 

where Z(/3 + \k) > 0 and 

7t(/3 + |/c) = —4 n 3 k 2 sn (/3 + \k) cn (/3 + \k) dn (/3 + \k) < 0 , (4.8) 

correspond to the energy gap (finite, upper forbidden band) E\ < E(a ) < E 2 , see Eq. (13.81) . in 
the spectrum of the Lame system. They are represented by real functions, and can be used in 
the Darboux transformation (12.81) . These functions, however, have infinite number of zeroes on 
the real line, and their use would produce a singular transformation. To resolve this problem, 
in the simplest case we can realize a second order (n = 2) Darboux-Crum transformation (12.101) . 
U 2 fl(x,t) = uo,o(x) — 2(log W (4>+(l), ^-(2)))xa:, by taking $+(1) = Cf $(x, t, ) + ^r$(x, t, ) 
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Figure 3: A pair of pulse solitons in the asymptotically periodic background of the KdV solution. 
Here /r = 1, k = 0.9, = 0.8, = 0.6, Cf = C 2 = 1. As in the case of asymptotically free 

background shown in Fig. [TJ a more deep defect propagates with a higher speed. Each of the 
two defects here supports a bound state in the lower forbidden semi-infinite band of the perturbed 
one-gap Lame system. 


and <£-(2) = C 2 ^{x, t, /3, + ) — -^4>(x, t, —with 0 < /3+ < P 2 < K and C+ 2 > 0. The non- 

singular potential U 2 ,o{x,t) supports two bound states in the spectrum of the Schrodinger operator 
L with energy values given by (13.5|) . which are in the gap, between the two continuous bands of 
the Lame system. The solution U 2 ,o(x,t) of the KdV equation, depicted in Figure 01 describes the 
propagation of the two solitons of the compression modulation type in the asymptotically periodic 
background. They are similar to the so-called grey solitons, see [5j [29i [30]. The described Darboux- 
Crum procedure can be generalized to obtain the potential supporting 21 bound states in the gap, 

U 2 i,o(x,t) = u 0 ,o(x) - 2(logfF($+(l),<f>_(2),...,<f> + (2/ - l),<f>-(2 l))) xx , (4.9) 

4> ( _ iy+ i (j) = C+$(x, t, pt) + (-1 y +1 ^(x, t, -/?+), j = 1,..., 21 , (4.10) 

0 < Pt < Pi < - < Pi < K, 0 < C+ < 00 . (4.11) 

The solution U 2 i t o(x,t ) has 21 compression modulation type solitons which move to the left in the 
stationary crystalline background. 

In definition of the wave functions (14. 5p and (14.101) . which are linear combinations of the 
eigenstates of the Lax operator P with opposite eigenvalues Tr(a) and — ir(a), the dependence 
on the evolution parameter t can be transferred from the common eigenstates of L and P into 
corresponding coefficients CJ and Cj~. In this way, function (14.101) can be presented equiva¬ 
lently in the form $(_ip+i(j) = CP(t)$(x,t = 0,Pf) + (~l) j+1 -^-^®(x,t = 0,-/3+), where 

Cj~(t) = Cf exp(7r(/3^")t), and 7r(/3+) is given by Eq. (14.81) with j3 + changed for /3~f. For (14.51) we 
have analogous equivalent representation with n(/3~ + /K / ) given by Eq. (14.31) . 
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Figure 4: A pair of solitons of the compression modulation type propagating to the left in a 
stationary crystalline background. Here (jl = 1, k = 0.9, /3± = 1, = 1-3, Cf = = 1. Due to 

the presence of the valence band in the spectrum of Lame system Lo,o, defects of this type have no 
direct analogs in the asymptotically free, homogeneous background case. 


To obtain the solution with odd number of solitons of the compression modulation type, there 
are the following possibilities. We can choose one of the arbitrary constants C' ? + and apply to the 
solution (14.91) the limit C+ —»• 0 (or c t —> oo). Such a limit results in sending the j -th soliton to 
x = —oo (or x = Too) without affecting the rest of solitons except inducing a global phase shift in 

(3+ f3+ 

them and in the background (equal tox->r —for Cj~ 0 and x —> x + for Cj~ —> oo) due 
to the nonlinear interaction in the KdV equation, and additional change in constants Cy, j' ^ j, 
depending on fJp and /3) 1 ", see [25] , Another option to generate the solution with odd number of 
solitons of the compression modulation type is to take /3f = 0 (in this case, < h+(l) oc sn(/rx|/c)), 
or @21 = K (then, <h_(2£) oc cn(/j,x\k)). One of the indicated states corresponds to the edge of 
the conduction band of the one-gap Lame system = 0), while another one corresponds to the 
upper edge of the valence band (/3^ = K). The same effect can be obtained just by applying the 
limit /3f —>• 0, or — > K in (14.91) . that results in eliminating the corresponding bound state from 
the spectrum of the Schrodinger system L = T U 2 i t o{x,t), see ref. [T2]. The issue of velocities 
for the described solutions will be discussed below, in Section m 

4.3 Mixed case 

In the mixed case, the Darboux-Crum transformation takes the form 

U 2 i, m (x,t) = u 0 , o(x) - 2(logW(4> + (l),4>_(2), ...,4>_(2/), ^+(1), J 7 _(2),..., {j))) X x ■ (4.12) 

This solution has 21 solitons of the compression modulation type, which move to the left, and m 
solitons of the potential well type, which propagate to the right. In the associated Schrodinger 
system L 2 z jm = — T U 2 i t m(x,t), the potential well and compression modulation type solitons 
support the bound states in the lower forbidden band and in the gap of the energy spectrum, 
respectively. In FigureO one can see how solitons of different types propagate in opposite directions 
over the stationary background. 
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Figure 5: The KdV solution with one pulse soliton (given by /3-f and shown in red) propagating to 
the right, and two solitons of the compression modulation type (characterized by /?(*" and /3^ and 
depicted in orange and green, respectively) moving to the left in a stationary crystalline background. 
More rapid compression soliton (/!+) supports a bound state of higher energy in the spectrum of 
the associated Schrodinger system L 2t i = — + u 2) i(x,t). At the moment t = 0, solitons are in the 
zone of a strong interaction, shown in violet, and are not well separated. The colour highlighting 
used here is rather conditional since the defects are transformed into a background asymptotically. 
The parameters are chosen so that the solution at t = 0 is symmetric with respect to the point 
x = 0. Here, // = !,& = 0.9, /Jj - = 0.9, = 1, ^ = 1-3, Cf = Cf 2 = 1. 


Here, the case with the odd number of the solitons of the compression modulation type can be 
obtained from (14.121) in the way described in the preceding Section 14.21 

4.4 On velocities and amplitudes of solitons in crystalline background 

For asymptotically free solutions discussed in Section [2~4l the multi-soliton solution can be presented 
in the form (I2.22jl . which is similar to the form of the second term on the right in ( 14 . 41 ) . ( 14 . 91 ) . or 
()4. 12 [) . Hyperbolic cosh and sinh functions in (I2.22|) are just the linear combinations of exponents 
exp(ttj(x — 4tt|t)) and exp(— Kj(x — 4tt|t)) with constant coefficients. As a result, when solitons are 
well separated (this happens for sufficiently large values of |t|), the solution (12.221) reduces to the sum 
of one-soliton solutions, u n (x,t ) ~ Y^j =1 — 2K 2 sech 2 (^Kj{x — xo j — 4K 2 t)^, with obviously identified 
velocities and amplitudes. Since the exponents in functions ( 14 . 11 ) and ( 14 . 71 ) are odd functions of the 
parameters /3~ and (3 + , respectively, functions ( 14 . 51 ) and ( 14 . 101 ) similarly to the hyperbolic functions 
are linear combinations of the exponents with arguments of the form (p(x — Vt) and —<p(x — Vf), 
but now with the coefficients which are certain periodic functions in x. As a consequence, unlike 
the case of KdV solutions over the asymptotically free background, it is not so obvious how to 
define amplitudes and velocities of solitons propagating in the crystalline background. Even when 
defects are well separated, their velocities and amplitudes are varied at each instant of time due 
to their nonlinear, position-dependent interaction with the oscillating background. Nevertheless, 
by analogy with the asymptotically free case, one can observe that the quantities /r 2 z 2 (/3 _ |fc) and 
/i 2 Z 2 (/3 + |fc) give us a relevant information on amplitudes of the well separated pulse and modulation 
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types defects, respectively, provided that for a well separated soliton defect the corresponding pre¬ 
exponential periodic factor in ((4.111 or (14.71) . associated with the crystalline background, will be 
less significant than the exponential factors there. This observation means a monotonic increase of 
the amplitude of the pulse type soliton with increasing energy modulus of the trapped by it bound 
state in the spectrum of the associated perturbed one-gap Lame system. Such a monotonic increase 
will be valid, however, starting from some sufficiently low energy value inside the lower forbidden 
band of the spectrum. For soliton defects of the compression modulation type the picture is more 
complicated since in the interval 0 < f3 + < K the function Z 2 ((3 + \k) takes zero value when j3 + tends 
to zero, then monotonically increases till some maximum value at some /3+ inside the indicated 
interval, and then monotonically decreases taking zero value at the edge /3 + = K. So, we can only 
conclude that the amplitude of the compression modulation type defect will tend to zero when 
the energy of the bound state trapped by it approximates the edges of the gap in the spectrum 
of the perturbed one-gap Lame system, and will take some maximal value at some energy of the 
corresponding trapped bound state inside the gap. With the same reservations one can say that 
the width of the defects is proportional to l//rz(/3 _ ) and l//iZ(/3 + ) for the pulse and modulation 
type solitons, respectively. 

To define the velocities of the solitons, we consider the case of a well separated defect of the 
pulse type, while for the compression modulation type defect the reasoning will be similar. 

The well separated pulse soliton with index j will have the form described by the one-soliton 
solution «o,i given by Eq. (14.41) with some parameter C{" and (3^ = j3j, i.e. this soliton will 
depend only on eigenstate J~+(j) with certain parameter Cj. Then we can identify the velocity 
of this soliton from the condition uo t i(x,t) = uo t i(x + A x ,t + A t ) to be valid for all x and t for 
which soliton remains well separated from other defects, where A x and A* are some constants. 
This condition has a solution with A x = nT, where T = corresponds to the period of the 
asymptotically periodic crystalline background and A t = n2K.z(/3~\k)/(fj.n((3~ + zK r |A:)), while n 
is any integer not violating the condition of separation of the soliton. From here we identify the 
velocity of the pulse soliton as V(/3“) = A x /A t , and obtain 


m) 


, 2 cn(/3 i \k)dn(/3 i |fc) 
M sn 3 (/3~| fc) 


(4.13) 


This corresponds to the velocity with which a fixed value of the argument of the exponent in (|4.1I) 
propagates in space and time, <p(x — Vt) = const. 

In a similar way, for the compression modulation type defect we find 




2 7„2 


-4 irk 


sn(/?+ \k)cn(0f |fc)dn(/3+ | k) 

Z(/3, + |fc) 


< 0 . 


(4.14) 


For f3j —> 0, the energy of the bound state trapped by the pulse defect tends to minus infinity, 
and the amplitude and velocity (14.131) of the soliton tend to infinity, while its width tends to zero. 
For j3j —>• K, the amplitude of the soliton tends to zero, its width tends to infinity, while the limit 
of (I4.13P is finite and equals 

lim V(/T~) = 4/r 2 £;' 2 — . (4-15) 

/3"->K J E 

Here E is the complete elliptic integral of the second kind, and as a consequence of the inequality 
relations > i m, one finds that the limit value (j4.151) is inside the interval (4/r 2 fc /2 ,4^). 
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Similarly, for the compression modulation type defect, for the limits when its amplitude tends 
to zero, its velocity tends to nonzero limits, 


lim V(j3f) 


-4/r 2 fc 2 


K 

K- E 


< 0 , 


lim V(/3+) 

/3+^K 


-4 n 2 k 2 k' 2 


K 

E - WK 


< 0 , 


(4.16) 


while its width tends to infinity. 


4.5 Galilean symmetry 

The described solutions can be modified by employing the Galilean symmetry of the KdV equation. 
This symmetry means that if u(x,t ) is a solution of the KdV equation ut = Guu x — u xxx , then 
uc(x,t ) = u{x — 6 Gt, t) — G is also a solution of the KdV equation for any choice of a real constant 
G. In the described solutions, the crystalline background over which the soliton defects propagate 
deforming it due to a non-linear interaction, was static. The application of Galilean transformations 
to the solutions will boost the defects and background and also shift vertically the solutions by the 
additive constant —G. This Galilean transformation will not change, however, the relative velocities 
between defects and their velocities with respect to the boosted crystalline background. As we shall 
see below, the Galilean symmetry of the KdV equation will play a crucial role in the construction 
of the solutions for the mKdV equation by means of the Miura-Darboux-Crum transformations. 


5 Unification of the KdV and mKdV equations by Miura- 
Darboux-Crum transformations, and exotic supersymmetry 

5.1 Miura-Darboux-Crum transformations 

The defocusing mKdV equation and the KdV equation, 

Vt Gv v x T v xxx — 0, Ut 6 uu x V u xxx — 0, (5.1) 

are related by the Miura transformation m- Namely, the substitution of = v 2 ± v x into the 
KdV equation gives 


uf - 6+ u± xx = (2v ± d x )(v t - Gv 2 v x + v xxx ). (5.2) 

The two relations in (|5.2jl mean that if v is a solution of the mKdV equation, then both vr 2 satisfy 
the KdV equation. On the other hand, if both u + and u~ are solutions of the KdV equation, and 
there exists v such that u± = v 2 ± v x , then v obeys the mKdV equation [25]. Unlike the KdV case, 
if v(x,t) is the mKdV solution, then —v(x,t) is also solution. 

From the chains of the Darboux-Crum transformations we know that the function V m (x,t) in 
(12.171) and the corresponding solutions of the KdV equation u m (x,t ) and u m -i(x,t) are related by 

V m (xi t) T {Vm(x, ty) x — U m —\ (x,t) X m , V m (x, t) ( Vm(x , t))® — u m (x, t) X m . (5.3) 

With taking into account the Galilean symmetry of the KdV equation, we displace x —>• x — 6A m t, 
and make a change u(x, t ) —» u(x — 6At, t) — A = u x (x , t). Eq. (15.31) will transform then into 

Vm(x,t) + (v m (x,t)) x = Idx,*), V^x^t) - (v m (x,t)) x = u^(x,t), (5.4) 

where v m (x,t) = V m (x — 6X m t,t). The v m (x,t) are therefore solutions of the mKdV equation. 


14 






Utilizing this observation and solutions of the KdV equation obtained by means of the Darboux- 
Crum transformations, we can construct infinite number of solutions for the mKdV equation. The 
KdV solutions from Section T2.31 will provide us with the mKdV solitons taking constant asymptotic 
values at infinity. On the basis of the KdV multi-soliton defects from the preceding Section, we will 
obtain solutions for the mKdV equation in the form of soliton defects propagating in the crystalline 
background. 


5.2 Exotic N = 4 nonlinear supersymmetry 

Before we proceed to the discussion of the solutions to the mKdV equation and their peculiarities, 
we note that on the basis of relations (15.41) . the corresponding solutions of the KdV and mKdV 
equations can be related by exotic N = 4 nonlinear supersymmetry. Besides an ordinary N = 2 
supersymmetry, it incorporates two non-trivial, Lax-Novikov integrals for the two Schrodinger 
subsystems associated with the KdV equation, as well as two additional supercharges to be higher 
order matrix differential operators. 

To reveal the exotic supersymmetric structure, we first note that an ordinary N = 2 supersym¬ 
metry in the form of superalgebra 


[C-miSm,a\ — 0, {Sm,a, Sm^b} — 2 C r 


, 6 = 1 , 2 , 


(5.5) 


with (J 3 identified as a Z 2 grading operator, is generated by the extended, matrix Schrodinger 
Hamiltonian C m and the associated Dirac Hamiltonian T> m , 


( ° \ v = ( 0 A m(m) \ 

V 0 L m {m) J ’ m V AnM 0 ) 


(5.6) 


where <S mj i = T> rn and S mi 2 = io:{D rn . The Schrodinger operators L m _\(rn) and L m (m ), L n (m ) = 
— -jA +u^ m (x, t), n = m — 1 , m, are related here by the Darboux (Miura) transformation generated 
by the first order operators A m (m) and dL(m), where A m {n ) = 37 — V m (x — 6A n t,t), and the 
solution of the mKdV equation v m (x,t) = V m (x — 6A m t,t) can be considered as a superpotential. 

Let us remind that for the quantum mechanical operators C m and T> m , parameter t corresponds 
not to the Schrodinger or Dirac evolution, but, instead, it is associated here with the coherent pe¬ 
culiar isospectral deformations of their potentials governed by the KdV and mKdV equations. Any 
Schrodinger system with multi-soliton or finite-gap potential can be characterized by a nontrivial 
integral of motion in the form of the Lax-Novikov higher derivative differential operator of odd order 
[21 EH EH- As a consequence, for the extended, matrix system C m composed from a pair of such 
peculiar Schrodinger subsystems related by a Darboux transform, the quantum mechanical N = 2 
supersymmetry associated with fermionic integrals S mA is extended up to the exotic N = 4 nonlin¬ 
ear supersymmetry incorporating two additional supercharges Q m ,a- Additional supercharges are 
composed from higher (even) order differential operators which intertwine the diagonal elements in 
C m , and together with building blocks of supercharges 5 mja (being the first order differential opera¬ 
tors A m (m) and A^ m {m)) factorize effectively the Lax-Novikov integrals of the extended Schrodinger 
system C m . For the discussion of a general structure of the exotic N = 4 nonlinear supersymmetry 
associated with finite-gap and soliton systems see refs. [isiMGsiGEiiira and further references 
therein. 

Specifically, here the extension of N = 2 supersymmetry up to the exotic N = 4 nonlinear 
supersymmetry happens as follows. We have started with the stationary solution ( 13 . 11 ) for the KdV 
equation to construct solutions u m (x,t). For the initial Schrodinger operator Lq = — + uo(x), 

uq(x) = uq,o(x), we can construct the first order operator ( 12 . 121 ) . A\(x,t) = 777: — v\(x,t), which 
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provides us with the intertwining relation A\Lq = L\A\, where L\ = — + u\(x,t). The Lax 
operator Pq = P{uq ) constructed on the basis of the stationary solution uq{x) satisfies equation 
(13.21) . [Pq,Lq] = 0. Being the third order differential operator, this is the Lax-Novikov integral for 
the Schrodinger system Lq. The intertwining relation A±Lq = L\A\ and the conjugate relation 
show that for the Schrodinger system L\, the fifth order differential operator P\(x,t ) = AiPoAj 
is the Lax-Novikov integral of motion, [P\,Li\ = 0. Applying then Galilean transformation with 
the parameter G = Ai to the KdV solutions uq(x) and ui(x,t), we obtain the fifth order operators 
Pi(l) = A\(1)Pq(x — 6 Aif)A.{(l) and Po(l) = Lq(X)Pq{x — 6 AH) = A{(l)Ai(l)Po(x — 6 Ait), which 
are integrals for Li(l) and Lo(l)> respectively. Being the product of the integral Pq(x — 6 AH) with 
the Schrodinger Hamiltonian Lq{\), in this case the integral Po(l) is reducible. We take it, however, 
to construct two bosonic integrals for the extended system C\, 

-Pl.l = ( P ° 0 (1) pi 0 (1) ) , Pi ,2 = <73^,1 , (5.7) 

with the matrix elements to be differential operators of the same order. 

To get the analogs of integrals (15.71) for a general case corresponding to the extended Schrodinger 
system described by C m , we have to change the fifth order differential operator A\PqA\ f° r 
A m P 0 AJ m that is a differential operator of order 2 m + 3 commuting with the Schrodinger oper¬ 
ator L m = —jfA + u m (x,t). Then we realize a usual Galilean transformation with the velocity 
—6A m to obtain the operator P m (m ) = A m (m)Po(x — 6A m t)Am(m), that is the Lax-Novikov inte¬ 
gral for the quantum system L m (m). In a similar way, as analog of L 0 (l)P 0 (x — 6X\t) we get 
operator = L m _i(rn)A m _i(m)Po(a; — 6 A m t)A| ra _ 1 (m) that commutes with P m _i(m), 

where A n (m) = A n (m)A n -i(m)...Ai(m). In this manner we obtain a pair of matrix operators 
Pm,i = diag (P m _i(m), P m (m)) and V m ,2 = GzP m ,i, which are the integrals for the extended 
Schrodinger system C m , and in addition to (15.51) . we get the relations 

[Pm,ai £m] = [Pm,a,'Pm,b] = [Pm,li<Sm,b\ = 0 . (5-8) 

At the same time, the commutator of P m ,2 with S mA , a = 1 , 2 , will supply us with a pair of new 
fermionic integrals Q m ,ai which are differential operators of even order 2 m, that also commute 
with P m , i, [P m ,i, Qm,b\ = 0 . So, the operator P rn ,\- like C m , is the bosonic central charge of 
the superalgebra. Together, two bosonic integrals P m ,a , a = 1 , 2 , allow us to distinguish the 
eigenstates corresponding to the four-fold degenerate eigenvalues inside the continuous allowed 
bands of C m . Also, they detect all the edge-states and all the bound states in the spectrum of C m 
by annihilating them. The fermionic integrals and Q, m ,a generate transformations between the 
‘up’ and ‘down’ eigenstates of the same eigenvalues in the spectrum of C m , and as usually, complex 
linear combinations of i and S m> 2 , and of Q m ,.\ and Q m ,2 will be creation and annihilation type 
operators for those eigenstates. The bosonic integral P m 2 generates a kind of rotation between 
supercharges S a and Q a mmm- 

The four fermionic supercharges S m ^ a and Q rn ,a and two bosonic integrals P m ,a together with 
the matrix Schrodinger Hamiltonian C m generate the exotic nonlinear N = 4 supersymmetry, 
whose super algebraic relations will contain the coefficients to be polynomials in the central charge 
C m . Such unusual nonlinear extension of supersymmetric structure related to integrable systems 
was discussed in different aspects in refs. [T2l [HI HU [28l [32j [33]. The anti-commutation relations 
for S mja in (15.51) . and similar relations for Q m ,a with right hand side to be a certain polynomial 
of order m in C m , together reflect the fact that the square of the Lax-Novikov integrals Pm,a , in 
correspondence with the Burchnall-Chaundy theorem GU, is a certain polynomial of odd order 
(2m + 1) in C m . 
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There also exist the cases of the systems, the explicit examples of which will be considered 
below, when the described structure of exotic supersymmetry can be reduced in the order of 
differential operators corresponding to the set of integrals V a and Q a . This happens when the 
Schrodinger potentials u^ n (x,t) and u^X 1 (x,t) are completely isospectral, and the ordinary N = 2 
supersymmetry generated by the first order supercharges S a in accordance with (15.51) is sponta¬ 
neously broken. Specifically, in the case when Schrodinger potentials u^ l n (x,t) and u'^_ l {x,t) 
have a difference in one bound state in the spectra of the systems L m {m) and L m _i(m), the 
spectrum of the corresponding Dirac Hamiltonian with the scalar potential v m (x,t) will contain 
one kink as a defect. If this is the case, there is no reduction in the structure of the exotic 
supersymmetry generators; the envelope of the corresponding oscillating eigenfunction ^ m (x) of 
the initial one-gap Lame system used in the Darboux-Crum construction will exponentially in¬ 
crease in both positive and negative directions of x. In contrast, when the envelope of 'L m (x) 
increases exponentially in one direction while in other direction exponentially decreases, the corre¬ 
sponding first order intertwining operator A m (m ) = X m _\ (m) = 37 — v m -i(x,t,X m ) generates a 
nonlinear shift in the already present soliton defects as well as in the background, without adding 
a bound state into the spectrum of L m (m) = in comparison with that of L m _i(m), 

X m -i(m)L rn -i(m ) = L m -i(m)X m -i(m). The superpotential u m _i(x,t,A m ) relating such a pair 
of isospectral Schrodinger systems can always be obtained from the appropriate superpotential 
Vm S \x,t ) corresponding to the associated irreducible extended Schrodinger system by one of the 
limits of the form lim Vm^ (x, t ) = ±u m _i(x, t, X m ) [XT, 125] . The spectrum of the corresponding 

C m —>0,oo 

Dirac Hamiltonian V m (supercharge S TOj i), with the scalar potential v m _i(x, t, X m ) is symmetric 
with a central gap between bound states or continuous bands, and it can contain defects of the 
kink-antikink type only, but never kink type defects. For such extended Schrodinger systems, the 
differential order of integrals V a and Q a will reduce by two. From the point of view of the limit 
of the associated appropriate irreducible system £m S \ we have —> C m , and then one can find 
that 'Pm S ' > — > C m V m , and a similar relation for supercharges Q a . The indicated reducibility of 
the matrix differential operators V a and Q a does not affect the nature of the supersymmetry: the 
corresponding extended Schrodinger system presented by the 2x2 diagonal matrix Hamiltonian 
operator C m is characterized by the exotic nonlinear N = 4 supersymmetry generated by two first 
order supercharges S a and by two fermionic integrals Q a being higher (even) order differential op¬ 
erators, and by two bosonic integrals composed from the Lax-Novikov integrals of the completely 
isospectral Schrodinger subsystems. 

6 Soliton solutions for the mKdV equation 

In the next two subsections, we discuss briefly the mKdV solutions corresponding to the multi- 
kink-antikink solitons propagating over the asymptotically free kink or kink-antikink backgrounds, 
and then we consider much more rich case corresponding to solutions in the crystalline kink and 
kink-antikink backgrounds. 

6.1 Multi-kink-antikink solutions over a kink background 

The employment of the multi-soliton solutions of the KdV equation constructed on the basis of 
the initial trivial solution u = 0 allows to find topologically nontrivial solutions v K (x,t) for the 
mKdV equation with the asymptotic behaviour ^(+ 00 ,f) = —v K (—oo,t) = const 7 ^ 0. Index 
K reflects here the kink-type nature of the solutions. To obtain the mKdV solutions with the 
indicated asymptotic behaviour, we should take u m with one more bound state in the spectrum 
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of the associated Shrodinger system in comparison with the (m — 1 ) bound states supported by 
the potential u m - In this case, the solution to the rnKdV equation will have a form of the 
multi-kink-antikink defect propagating over the kink. The kink-antikink perturbations will always 
have amplitudes smaller than the kink amplitude, characterized by the parameter K m , due to the 
ordering 0 < < ... < n m . The kink-antikinks ( j = 1 ,... ,m — 1 ) and kink (j = m) propagate 

to the left with the velocities Vj = 4 k 2 — 6 k^. So, the kink’s speed is the lowest, being equal to 
2 k^, while the kink-antikink solitons with lower amplitudes have higher speeds. Such solutions 
generalize the mKdV kink solution = Ktanh/c(x + 2n 2 t), and analytically they are given 

by 

Vm-i(x,t) = (x + 6k^M), (6.1) 

where, in correspondence with (12. 220 . 


V*(x,t) = ftm-iOM) a m = - (log W (cosh X 1 ,sinh X 2 ,...,f(X m ))) , ( 6 . 2 ) 


/(X“) = sinh X~ if m is even, /(X“) = coshA“ if m is odd, and X~ is defined in (12.22p . 
Functions ±u^(x,f), representing solutions for the defocussing mKdV equation, are of the kink 
type with m solitons that deform in their propagation the moving kink (or antikink) background 
without overpassing its asymptotes. An example of such solution is represented in Fig. [ 6 l 


6.2 Multi-kink-antikink solutions over the topologically trivial background 

The multi-kink-antikink solutions propagating over the topologically trivial background with 
asymptotic behaviour v(—oo,t) = v(+oo,t) = const / 0 can be obtained from the KdV solu¬ 
tions u m and u m -\ which, as the Schrodinger potentials, support the same number [(m — 1), 
m = 2,...] of the bound states with coinciding energies. In such a pair, potentials u m and u m -\ 
are related by the Darboux transformation that displaces solitons | T2| [28] . They generalize the 
simplest kink-antikink solution for the mKdV equation, 

v 0 (x,t) = -K 1 , (6.3) 

which corresponds to the nonzero mass term rn = n\ of the free Dirac Hamiltonian operator Dq, 
see (15.61) . 

The simplest generalization of (16.31) is given by 

v\{x , t ) = «q tanh (ki(x + vt)) — n\ tanh (ki(x + vt) — A) — n\ coth A , (6.4) 

where u = 6^2 ~ 4/cf > 0, A = ^ log > 0, and K 2 > > 0. The mKdV solution (16.41) is 

related to the mutually displaced one-soliton KdV solutions u\(x, t ) and U 2 (x, t ) via the Miura trans¬ 
formation, u\ = + = k?; — 2 /vj ! sech 2 (k±(x + id)), U 2 = v\ — v[ = — 2n\sec\\ 2 (ki(x + vt) — A). 

Analytic form of this type of the solutions is given by relations of the same form (16.ip . (| 6 . 2 I) but with 
f(X ~) changed here for f(X~) = exp A“. Such solutions can be obtained from the topologically 
nontrivial solutions (| 6 . 2 I) by taking there a limit xo m +oo or —oo, that corresponds to sending 
the soliton (kink) associated with the non-degenerate zero energy in the spectrum of T> to +oo or 
—oo. Thus, the velocities of the remaining kink-antikink solitons with j = 1,... ,m — 1 are the 
same, Vj = 4k 2 — 6 k^, as in the solutions with the kink, indexed there by j = m, and its traces are 
still present here by restricting the amplitudes of the kink-antikinks and their speeds. An example 
of such type of solutions is shown in Fig. [71 
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Figure 6 : Solution of the rnKdV equation representing the kink-antikink (characterized by the k\ 
parameter) propagation over the kink (given by K 2 ) background, both moving to the left. When 
the soliton, which is faster, is to the right from the kink center (see the t = —3 instant), it has 
a form of the antikink-kink perturbation. Overtaking the kink center, it flips and transforms into 
the kink-antikink perturbation (t = 3 and t = 6 ). Here ki = 1, «2 = 1-1, xqi = 0, i = 1,2. 
The change v —» — v provides a solution over the antikink background with a mutual transform 
of the antikink-kink and kink-antikink perturbations. The kink has a bigger amplitude 2^2 equal 
to the distance between asymptotes, and this value corresponds also to the size of the central gap 
in the spectrum of the associated Dirac Hamiltonian operator shown on the right. Filled circles 
correspond to non-degenerate energy levels at the edges of the doubly degenerate continuous parts 
of the spectrum, while blue squares indicate non-degenerate bound states inside the spectral gap. 
Zero energy bound state is associated with the kink, the two other bounds states are associated 
with the kink-antikink perturbation in the rnKdV solution. 


6.3 Kink-antikink pulse type defects in a crystalline kink background 

On the basis of the stationary cnoidal solution for the KdV equation from Section [3l one can 
construct diverse types of solutions for the rnKdV equation, some of which are topologically trivial, 
while others have a nontrivial topological nature. If we use as u m -\ a solution u m -i t 0 of the KdV 
equation that contains only the potential well soliton defects in the periodic background, and take 
u m = Um ,0 of the same type but with one pulse soliton defect more, then the associated solution v m 
of the mKdV equation will describe multiple kink-antikink perturbations propagating in a moving 
crystal kink, see Figure [H The kink here is the defect of the greatest amplitude. Together with 
other soliton defects, it propagates to the left, see Section 16.101 below. In this case the solutions 
can be presented in the form 

v o,m-1 (®> t) = V 0 K m (x - 6A(/3“ + iK')t, t) , (6.5) 

where 

Vo K m(- M) = = - (log W(F + (1) , (2),..., J r ( _ 1) m +1 (m))) x , 

( 6 . 6 ) 
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Figure 7: Solution of the mKdV equation corresponding to two kink-antikink perturbations propa¬ 
gating to the left. Here K\ = 1, K 2 = 1.2, = 1.2+10 , xo , = 0, i = 1,2, 3. 2^3 corresponds to the 

size of the gap in the spectrum of the Dirac Hamiltonian operator, K 2 defines a bigger kink-antikink 
soliton; the closeness of its value to K 3 prescribes this soliton to be higher and wider, and defines 
the energies ±a/k 2 — K 2 °f the two bound states supported by this defect. k\ defines parameters of 
the smaller kink-antikink soliton and the energies ± \J /c 2 — k 2 of the two bound states supported by 
it. In contrast with the KdV solutions with a free background, the defects with smaller amplitude 
propagate with higher speed. A nonzero asymptotic value of this type solution corresponds to a 
constant mass term in a free massive Dirac Hamiltonian operator. 


and A (a) = /i 2 (dn 2 ( a\k ) — 1(1 + A/ 2 )). The spectrum of V is symmetric, with two finite, and two 
semi-infinite allowed bands. It contains a finite number of bound states in central gap, and one of 
these bound states is exactly in the center (E = 0), while two other gaps are unoccupied. 

6.4 Multi-kink-antikink pulse type defects in a kink-antikink crystal back¬ 
ground 

Let us take a solution u m _i j 0 of the KdV equation as a solution u m -\ in (15.31) . and choose u m in 
the form of the solution of the same type, u m = u m ~i t 0 , but displaced by means of the Darboux 
transformation. Then we get the associated solution v m of the rnKdV equation in the form of 
multi-kink-antikink defects propagating in a crystalline background. In this case, again, both the 
crystalline background as well as the pulse defects will propagate to the left. Such type of solutions 
can be presented in the analytical form 

vo,m-i = Vo,m(x - 6A(/3“ + iK')t, t ), (6.7) 

where 

Vo,m(x,t) = (log W (J 7 4 .( 1 ), j 7 -(2), ..., J 7 (m — l)))x 

- 1 ),F(x,t,/3~))) x . (6.8) 

Function F(x, is defined here by Eq. (|4.1I) . The spectrum of T> is symmetric, with two finite 

and two semi-infinite allowed bands. It has finite even number of the bound states in the central 
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Figure 8 : A kink-antikink (given by (3^) and a kink (given by j3^) type defects propagating over 
a kink-antikink crystalline background. This mKdV solution with n = 1, k = 0.999, (3= 1.7, 
(3 2 = 1-5 and C~ = 1, i = 1,2, is a somewhat analogous to that presented in Figure El The kink 
perturbation and the kink-antikink in the form of pulse (pul) soliton defect as well as the kink- 
antikink crystalline background (bg) move to the left, and magnitudes of their velocities (speeds) 
are subject here to the relation 0 < \Vkink\ < \Vpui\ < | 1, see Section IB. 101 below: the pulse soliton 

defect flips when passes from one to another side of the crystalline kink defect. Notice that in 
contrast to the present case, for the mKdV solution depicted in Figure El the background is given 
there by constant asymptotes for which a state of motion is not defined. The size of a central 
gap in the spectrum of the associated Dirac Hamiltonian operator shown on the right is given by 
2/r|dn(/3rj" + zK'l/c)|. The energies of the two bound states trapped by the kink-antikink defect are 

given by E±(a) = ±/ry^dn 2 (a\k) — dn 2 (/3^~ +iK.'\k) with a = /3{~ + iK', while the zero energy 

value, E±(a = (3% + iK')) = 0, corresponds to a unique bound state trapped by the kink defect. 
The edges of the allowed bands correspond to a = 0, K, K + iK'. 


gap, and so, in contrast with the previous case, there is no bound state of discrete zero energy in 
the center of the gap. As in the class of solutions discussed in the previous Section, the symmetric 
non-central gaps are empty of bound states. The mKdV solutions of this type can be obtained from 
those described in the preceding Section by sending the kink to plus o minus infinity; the concrete 
analytic form (16.811 corresponds to taking the limit C~ —>• oo in the solution (16.5|) . This case is 
illustrated in Fig. El 

6.5 Multi-kink-antikink compression modulation defects over a kink in a kink- 
antikink crystal background 

There are no nonsingular solutions of the mKdV equation associated with the KdV solutions which 
in the Darboux-Crum construction use only the states from the gap. To get nonsingular solutions, 
one can employ in the construction of u m a state from the forbidden lower band of Lame system 
in addition to the states from the gap employed for the construction of or to generate a 

nonlinear displacement by means of the Darboux transformation. In the first case we obtain a 
solution for the mKdV equation which contains a kink in a background, while in the second case 
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Figure 9: A kink-antikink pulse (given by /3^) over a kink-antikink crystalline background. Here 
H = 1, k = 0.9999, = 1.5+10 -12 , /3^ = 1.5, and C~ = 1, i = 1, 2. When the values of and (5^ 

are closer, the kink-antikink defect in the background of the kink-antikink crystal is more notable. 
Both the kink-antikink pulse and the kink-antikink crystal background move to the left, and the 
speed of the latter is higher than that of the defect. The size of a central gap in the spectrum of 
the associated Dirac Hamiltonian operator shown on the right is given by 2/r|dn(/3^" +iK , |A:)|. The 

energies of the bound states trapped by the defect are E±(a) = ±/i-^dn 2 (a\k) — dn 2 (/3^~ + iK.'\k) 

with a = /3j~ + iK 7 , and they are represented by the blue rectangles. The non-degenerate energy 
values of the edges of the allowed bands correspond to a = 0, K, K + iKJ. 


there will be no such kink in the structure of the mKdV solution. The solution of the first indicated 
case takes the form 

v 2i,o = v u,\ ( x ~ 6A(£f + iK%t ), (6.9) 

where 


<1 (M) = (log VF(<3? + (1), <h_(2),..., < f>+(2/ — 1), <h__(2/), J 7 4.(l))) x 

-(log W($ + (l), < h-(2),..., <h+(2 1 - 1), <h_(2 l))) x • (6-10) 

In this solution, the defects of the compression modulation type propagate over a kink defect, 
which, in turn, propagates over the the crystalline background by deforming it. An example of such 
solution is shown in Fig. [TOJ Here, the kink, the background, and the compression modulation type 
defects move to the left. The shown solution possesses an even number of compression modulation 
defects. The solutions with odd number of modulation defects can be obtained by applying the 
procedures explained earlier for the KdV equation. The spectrum of 0 is symmetric, with two 
finite allowed bands and two semi-infinite allowed bands. Besides, it has finite number of bound 
states in the external gaps, and one bound state of zero energy in the central gap. 
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Figure 10: The mKdV solution with two kink-antikink modulations (given by (3+ and /3^) and a kink 
(given by /3j” and corresponding to the KdV pulse soliton) propagating in a kink-antikink crystal 
background. Here n = 1, k = 0.9, /3{~ = 1.8, /Sj 1 " = 1, /3^~ = 1.3 and Cj - = Cj 1 " = = 1, and the 

velocity magnitudes of the solitons of the compression modulations (mod) type, of the background 
(bg), and the kink soliton are subject to inequalities |V m0( z| > \Vbg\ > \Vkink\ > 0, see Section 16.101 
below. The size of the central gap in the spectrum of the associated Dirac Hamiltonian operator, 
shown on the right, is equal to 2/x|dn(/3j" + zK 7 17c)|. The energies of the bound states trapped 

by the kink-antikink modulations are given by E±(a ) = ±fi-^dn 2 (a\k) — dn 2 (/3j~ + iK.'\k ) with 

a = (3i , f3 2 , while the bound state of zero energy (a = /3{~ + ?'K') is trapped by the kink defect. 
The edges of the allowed bands in the spectrum correspond to a = 0, K, K + zK 7 . 


6.6 Multi-kink-antikink compression modulation solutions in a kink-antikink 
crystal background 

If the difference between u m and u m -i is a nonlinear displacement generated by a Darboux trans¬ 
formation, and if these solutions of the KdV equation are of the compression modulation type 
defects, then the associated solution of the mKdV equation will be a kink-antikink crystal back¬ 
ground propagating to the left, in which we have kink-antikink defects of the modulation type also 
moving to the left. Analytically such solutions are given by 

V 2 i,o = V 2 i,i(x - 6A(/3f + iK')t, t ), (6.11) 


where 


Vu,i(x, t ) = (log W (3>+(l), M2), - 1), M2 l),F(x, t, Pi))) x 

— (log W (<h + (l), <I>_(2),..., M2 1 ~ 1), ^-(2/)))* • (6-12) 

Notice a difference in the last argument of the Wronskian in the first term on the right in (16.121) in 
comparison with (16.101) . It reflects the fact that the solutions of the present type can be obtained 
from those discussed in the preceding Section by sending the kink to infinity, by taking the limit 
C7f — > oo in (16.10|) . 
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Figure 11: The mKdV solution with two kink-antikink modulations (given by /3f and /3^") prop¬ 
agating over a kink-antikink crystal background (described by the parameter Pi). Here /j = 1, 
k = 0.9, (3 = 1.8, (3f = 1, P 2 = 1 - 3 , C ±2 = 1 - At t = 0, the vertical axis (x = 0) is shifted 
a little bit to the left with respect to the symmetry axis of the graphic. This is due to the limit 
Cf 00 applied to the solution shown on Figure [TO] the present solution. This limit corresponds 
to sending the kink to infinity and generates the displacements (phase shifts) in the background 
and in the remaining defects. The size of the central gap in the spectrum of the associated Dirac 
Hamiltonian operator is equal to 2/t|dn(/df + iK'|fc)|. The energies of the bound states trapped by 

the kink-antikink modulation defects are given by E±(a) = ±fi-^dn 2 (a\k) — dn 2 (/Lf +iK.'\k) with 
a = (3f, @ 2 - The edges of the allowed bands correspond to a = 0, K, K + IK'. The compression 
defects are more rapid than the background. 


In Figure [TT] it is shown the case with even number of the compression modulations defects. 
The case with the odd number of such defects can be obtained by the procedures explained earlier 
for the KdV system. The spectrum of X> 2 i t 0 is symmetric, with a finite number of bound states 
appearing in non-central gaps, and with no bound states in the central gap. 

6.7 Multi kink-antikink modulation solitons in a kink crystal 

There is a special case which can be obtained as a limit from the solutions discussed in the preceding 
Section. It is generated by the choice of the solution u m -1 of the KdV equation which contains 
only compression modulation defects, while the additional state T oc dn(/j,x\k) employed for the 
construction of u m is at the edge of the lower forbidden band of the associated one-gap Lame 
system In this case the crystalline background in the mKdV solution is centered (vertically) 
in zero in contrast with the previously considered cases where it was displaced up or down. This 
centered crystalline background is known as the kink crystal solution that appears in the Gross- 
Neveu model OH H6]. In comparison with the previous cases, here in the spectrum of T> the central 
gap (together with bound states there) disappears and two finite continuous bands merge into one 
central allowed band centered at zero. So, in this case one can have defects only of the compression 

2 This is the ground state of the Lame system at the lower edge of its valence band. 
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modulation type, see Fig. H2J Here, as in the previous cases, the indicated soliton defects move 
to the left like the kink crystal propagating with the velocity 6A(K + iK 7 ), see Section f6. 10 1 below. 
Analytic form of such type mKdV solutions is given by 


V 2 i C = V 2 f c (x-6\(K + iK')t,t), 

where 


(6.13) 


v£ c (x,t) = (log W (d> + (l), $-(2),..., < h+(2/ — 1), <h_(21), dn (/j,x\k))) x 

— (log W (<h + (l), <h_(2),..., <&+(2l — 1), <£>^(2/)))^ . (6.14) 

In the simplest case l = 0 here the solution is just the kink crystal 


v£ c (x,t) v(p% c ) 



Figure 12: The rnKdV solution with two kink-antikink modulations (given by /3f and fi 2 ) over 
a kink crystal. Here, fi = 1, k = 0.9, /3f = 1, /3 2 = 1.3, /3j~ = K, = Cf 2 = 1- In this 

case there is no central gap in the spectrum of the associated Dirac Hamiltonian operator. The 
energies of the bound states trapped by the modulation defects are given by E±(a ) = ±/rdn(a|fc) 
with a = Pi , P 2 i w hde the values a = 0, K correspond to the edges of the allowed bands. In 
the kink crystal background, only the modulation type kink-antikink defects can exist, and they 
propagate to the left more rapidly than the kink crystal background. In this configuration, the 
velocity magnitude of the background, |Vb s | = 2/i 2 (l + A/ 2 ), is minimal in comparison with that in 
other types of the mKdV solutions. 

v o C (x,t) = (log dn (fi(x — Vb g t)\k)) x (6.15) 

propagating to the left with velocity Vb g = —2/r 2 (1 + A/ 2 ) < 0. 

6.8 Mixed multi-kink-antikink solitons over a kink in a kink-antikink crystal 
background 

In a more general case, one can have both types of defects, compression modulations as well 
as pulse solitons, propagating over a moving crystal background. In dependence on whether the 
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difference between u m -1 and u m solutions of the KdV equation is a pulse type soliton or a nonlinear 
displacement, there will appear or not a kink defect in the kink-antikink crystal background. In 
the case when the indicated difference is a pulse type defect, the solutions over the kink crystalline 
background take the form 

«a,m—1 = V £,m( x - 6A (/ 3 m + *) , (6.16) 

where 


v£ m (x,t) = (log W ($+(l), <h_(2),..., <I>„(2/), J 7 4 _( 1 ), J 7 _(2),..., J 7 ’(-ijm+i (m))) x 

— (log W (<h + (l), 3>_(2),..., (21), J 7 ^(1), J 7 _(2),..., J 7 (_ 1} m (m — l)))j; (6.17) 

is a modulated kink. In this solution, the kink, the pulse and modulation type defects, as well as 
the crystalline background move to the left. 

The shown in Fig. [13] solution possesses even number of compression modulation defects. 
Again, solutions with odd number of defects of this type can be obtained by means of any of the 
procedures discussed in Section 14.21 The spectrum of the Dirac Hamiltonian operator T> in this 
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Figure 13: The rnKdV solution with a kink-antikink pulse (given by and highlighted in red) and 
two kink-antikink modulations (given by (3f and /3^~, and shown in orange and green, respectively) 
propagating over a kink defect (/3^~), with all these defects propagating in a moving kink-antikink 
crystal background. Here n = 1, k = 0.999, / 3 ± = 1.7, j3^ = 1.5, (3^ = 1, = 1.3, and 

Ci 2 = 1, and the magnitudes of velocities of the defects and background are subject to inequalities 
|V mo d| > |Vbg| > |Fp u i| > |Vkink| > 0. The size of the central gap in the spectrum of the associated 
Dirac Hamiltonian operator is 2/r|dn(/3^~ +iK , |fc)|. The energies of the bound states trapped by the 

modulation defects are given by E±(a) = ±fi-^dn 2 (a\k) — dn 2 (/3^" + iK.'\k) with a = /?+, (3% ■ The 
energies in the central gap, which are marked by two blue squares symmetric with respect to zero 
level, are given by a = + zK 7 and correspond to the bound states trapped by the kink-antikink 

pulse; the zero energy value, E±(a = ^ +7K')) = 0, corresponds to a unique bound state trapped 
by the kink defect. The edges of the allowed bands correspond to a = 0, K, K + i K 7 . 

case is symmetric. Besides the continuous bands shown on the figure, it contains a finite number 
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of bound states in external and central gaps. Besides, there is an additional bound state of zero 
energy in the center of the central gap, which is associated with the crystalline kink. 

6.9 Mixed multi-kink-antikink defects in a kink-antikink crystal background 

If, unlike the previous case, the difference between u m -± and u m is given by a nonlinear displace¬ 
ment, there is no kink structure in the corresponding mKdV solution given by 

V2i,m-1 = v 2 i t m(x - 6A (/ 3 ~ + iK')t, t) , (6.18) 

where 

v 2 i,m(x,t) = (logW($+(l),$_(2), ...,$„(2/),^+(1), J 7 _(2), - 1 ),F(x,t,/3~))) x 

-(logW($ + (l),$_(2),...,$_(2Z),J-+(l),J-_(2),...,^ ( _ 1)m (m-l))) x (6.19) 

is a kink-antikink crystal background displaced from zero in vertical direction and perturbed by 
propagating in it mixed multi-kink-antikink defects. The example of such a solution is depicted 
in Fig. [TU The symmetric spectrum of T> has a finite number of bound states in the central and 
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Figure 14: The rnKdV solution with a kink-antikink pulse (given by and shown in red) and two 
kink-antikink modulations (given by f3+ and /3^~, and highlighted in orange and green, respectively) 
propagating to the left in a moving to the left kink-antikink crystal background. Here /i = 1, 
k = 0.999, /df = 1.5 + 10“ 10 , /3^ = 1.5, = 1, = 1.3, C± 2 = 1> and the velocity magnitudes 

are subject to inequalities |V mo d| > |Vb g | > |V pu i| > 0. The energies in the spectrum of the 
associated Dirac Hamiltonian operator are given by expressions similar to those in Figure EB 

external gaps, with no state of zero energy in the center of the central gap. The mKdV solutions of 
this type can be obtained from those discussed in the preceding Section by sending the kink defect 
to infinity. 

6.10 Velocities in the mKdV solutions with crystalline backgrounds 

The mKdV solutions v m -\ are obtained on the basis of the Darboux construction of KdV solutions. 
The necessary step of the procedure, as we have seen, involves the Galilean transformation, x —> 
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x — 6A m t. The boost parameter G = 6A m of this transformation is given by the energy A m that 
corresponds to the nonphysical eigenstate of the one-gap Lame system Lo,o which is used to obtain 
the potential u m with one more bound state in comparison with the partner potential u m —\ . Or, 
this X m corresponds to the eigenstate with the help of which the potential u m is obtained from the 
completely isospectral potential u m -± by means of the nonlinear Darboux displacement. In both 
cases A m = A(/3“+iK'), and /3“ in the first case is associated with the mKdV kink supporting bound 
state of non-degenerate discrete zero energy eigenvalue in the spectrum of the Dirac Hamiltonian 
operator V, or, if there is no discrete zero energy value in the spectrum of V, /3“ can be associated 
with the kink defect sent to infinity. This /3“ corresponds to the eigenstate of the minimal energy 
from the lower semi-infinite forbidden band in the spectrum of Lo,o which is used in the Darboux- 
Crurn transformations to construct u m , or to the eigenvalue of the ‘virtual’ state from that band 
associated with the ‘kink sent to infinity’. The crystalline background in the obtained rnKdV 
solution propagates as a result to the left and its velocity is 

V bg = 6A(/3" + iK') = —2fi 2 (3cs 2 (P~\k 2 ) + 1 + k' 2 ) < 0 . (6.20) 

The background propagates with a minimal speed in the class of the rnKdV solutions with the kink 
crystal background, which were discussed in Section [6.71 In that case = K corresponds to the 
edge-state dn(/j,x\k ) of To,o> and, therefore, Vb g = —2/i 2 (l + k /2 ) < 0. 

If the mKdV solution has the pulse type defects different from kink, their velocities are given 
by 

V pul (pr) = V(pr) + 6\(p-+iK')<0, K>pr>p~, j = ( 6 . 21 ) 

with V((3J) given by Eq. (Id.hip . Analogously, the velocities of the modulation type defects, if they 
are present in the mKdV solution, are 

Vmod(Pf) = V(P+) + 6\(P~+iK') < 0 , 0<P+<K, ( 6 . 22 ) 

where V(pP) are given by Eq. (14.1411 . In the case if the mKdV solution has a kink defect, its 
velocity is defined by relation of the form (I6.2ip but with j = m, i.e. 

Vkink = V(P~) + 6A(/3" + iK') < 0. (6.23) 

One can see that in the mKdV solution, the speeds of the background and defects, some types 
of which can be absent in the solution, are ordered according to the inequalities |V mo <i| > \Vbg\ > 
|kp U z| O 0. 


7 Discussion and outlook 

We have constructed solutions to the KdV equation with an arbitrary number of solitons in a 
stationary asymptotically periodic background. In this case there exist two types of solitons: 

• potential well defects (pulses), which propagate to the right, 

• compression modulation defects, which move to the left. 

These solutions asymptotically have a form of the one-gap Lame potential but subjected to the phase 
shifts x —>• x — x^ = ±-j(JV + Ylj Pj~) f° r x Too with respect to the stationary solution 
m- The asymmetry in propagation of the two types of soliton defects is valid in the case of a 
stationary background. If we apply Galilean transformations to the KdV solutions, we obtain new 
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solutions, for which, in general case, the described propagation asymmetry of the defects over the 
now moving background will be changed. However, this does not change the picture of the relative 
motion: pulse defects always will propagate to the right with respect to the asymptotically periodic 
background, while modulation type defects will move to the left with respect to non-stationary 
crystalline background. The interesting peculiarity we also have observed in the constructed KdV 
solutions is that in the limit cases when amplitudes of the pulse and compression modulation defects 
tend to zero, the limit values of the velocities of defects with respect to the crystalline background 
are nonzero. 

For the mKdV equation, we have constructed the following solutions from the obtained KdV 
solutions by means of the Miura-Darboux-Crum transformations: 

• solutions with a kink crystal background, in which there can exist only solitons in the form 
of the compression modulation type defects, 

• solutions with a kink-antikink crystal background, in which there can exist kink-antikinks in 
the form of the pulse and/or compression modulations type defects. 

In a kink-antikink crystal background there also can appear a topological defect in the form of the 
kink (or antikink, if we use a symmetry of the mKdV equation by changing v(x,t) for —v(x,t)) 
which always is related with the KdV pulse type defect that traps the bound state with the lowest 
energy in the lower forbidden band in the spectrum of the associated perturbed one-gap Lame 
system. 

Unlike the KdV, the mKdV equation has no Galilean symmetry, and the velocities of the defects 
and asymptotically periodic background in the solutions we constructed have an absolute character. 

In the mKdV solutions, all the defects and crystalline backgrounds move to the left, and the 
velocity magnitudes of the kink-antikink defects of the modulation (mod) and the pulse (pul) types, 
and the velocity magnitudes of the background (bg) and the kink are subject to the inequalities 
|Vmod| > |Vfe g | > \V pu i\ > |Vfcj n fc| > 0. Thus, with taking into account the sign of the velocities, 
we have for the defects and backgrounds in the mKdV solutions, similarly to those in the KdV 
solutions, ( Vp U i — Vb g ) > 0 and ( V mo d — Vbg) < 0. At the same time, the velocity magnitude (speed) 
of the kink defect, if it is present in the mKdV solution, always has a minimal value in comparison 
with other velocity magnitudes. 

The presence or absence of the kink in the mKdV solution is detected by N = 2 supersymmetry 
of the associated extended Schrodinger system. It is generated by the first order supercharge 
operators S a . a = 1,2, having a nature of the Dirac Hamiltonian operators with a scalar potential. 
When kink is present or absent, the supercharges possess or not a zero mode, and the N = 2 
supersymmetry is unbroken or broken. The supersymmetry detects also the case of the mKdV 
solutions with the kink crystal background. For such solutions, the kernel of supercharges is two- 
dimensional. Unlike the solutions with kink, the corresponding zero modes in this case are given 
by not normalizable but periodic states El- The position of the bound states in gaps in the spectra 
of the supercharges defines also the type of defects present in the corresponding mKdV solution. 

We have showed that the N = 2 supersymmetry constitutes a part of a more broad, N = 4 type 
exotic nonlinear supersymmetry, which includes in its structure two bosonic generators composed 
from the nontrivial, Lax-Novikov integrals of the pair of the Schrodinger subsystems. These bosonic 
generators are higher derivative differential operators, one of which is the central element of the 
superalgebra. Besides, exotic supersymmetric structure contains additional pair of supercharges 

3 These zero modes are constructed from Darboux-dressed edge states dn(fix\k) and dn {fix + K|fc) of the pair 
of mutually shifted in the half-period Lame systems by applying to them the Galilean boost with velocity Vbg = 
- 2^ 2 (1 + fc ' 2 ). 
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being matrix differential operators of the even order. The additional integrals appropriately reflect 
the peculiar nature of the extended Schrodinger system associated with the pair of the KdV solutions 
by detecting all the bound states and the band-edge states in the spectrum, as well as distinguish 
the eigenstates corresponding to the fourth-fold degenerate energy values inside the allowed bands. 

The both, KdV and mKdV equations are invariant under simultaneous inversion of t and x 
variables, but not under reflection of t (or, of x) only. The stationary trivial solution u = 0, 
on the basis of which we construct soliton solutions over the asymptotically free background, is 
invariant under separate inversions in t and x. The same is true for the stationary cnoidal solution 
in the form of the stationary one-gap Lame potential. So, the ‘initial data’ in the construction 
are invariant under t- inversions, and anisotropy of evolution as well as the chiral asymmetry of 
the non-stationary solutions is rooted in the anisotropy of the equations themselves. Namely, in 
our construction, though the ‘seed’ KdV solutions are time-inversion invariant, the solutions of the 
auxiliary problem in the Lax formulation for the KdV equation, which are generating elements 
for the Darboux-Crum transformations, break this symmetry. In the case of the asymptotically 
free background, in the indicated solutions of the auxiliary problem the dependence on time enters 
universally in the form of arguments x — Vjt with Vj = 4 k j >0. As a result, all the solitons in non¬ 
singular KdV solutions move to the right. For the soliton solutions over an asymptotically periodic 
background, the third order Lax operator of the initial Lame system has eigenvalues of different 
signs on the upper and lower horizontal borders of the fundamental ct-rectangle which correspond to 
the lower forbidden band and the gap. It is this sign asymmetry that finally is responsible for chiral 
asymmetry in propagation of the KdV solutions in the form of pulse and compression modulation 
defects over a crystalline background. The asymmetry of the rnKdV solutions is inherited from 
that for the KdV solutions. 

Since the KdV and mKdV solutions, and particularly those associated with the Lame quan¬ 
tum system, find many diverse applications in a variety of different areas of physics ranging from 
hydrodynamics, plasma physics, and optics to hadron physics and cosmology E 0 El IH El IS 
El El EH EH ESI EH So], it would be very interesting to find where the obtained new solutions 
could show up themselves. They could appear in the form of perturbations of different nature, 
which would propagate in a nonlinear media with different velocities and reveal chiral asymmetry 
in their dynamics. Another peculiarity which could be associated with the described solutions is 
the existence of nonzero bounds for the velocity of the defects with disappearing amplitudes. 

It would also be interesting to consider a generalization of the approach employed here by using 
Darboux transformations for the first order Hamiltonian operator of a (l+l)-dimensional Dirac 
system instead of the second order Schrodinger operator. In this way one could get finite-gap Dirac 
Hamiltonian operators of the form 


V = 


Ux) ~^ + V 1 (x)\ 
£ + Mx) ~V 2 (x) ) 


(7.1) 


with asymmetric spectrum. The corresponding stationary potentials V\ ^(x) could then be promoted 
to solutions in the form of the twisted kinks and twisted kink-antikinks |18[ [34 , 1411142] in a periodic 
background for nonlinear Schrodinger equation belonging to the Zakharov-Shabat - Ablowitz-Kaup- 
Newell-Segur hierarchy. 

In the supersymmetric quantum mechanical structure we discussed, the N = 4 refers to the 
number of supercharges appearing in the extended Schrodinger system. As the extended system is 
composed from a pair of the perturbed one-gap Lame systems, one could expect the appearance of 
only two supercharges as it happens in supersymmetric quantum mechanical systems of a general 
nature [43]. The peculiarity of the considered systems consists in their finite-gap nature, and it 
is this property that is behind the extension of the usual N = 2 supersymmetry up to the exotic 
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N = 4 nonlinear supersymmetric structure that incorporates the pair of Lax-Novikov integrals 
of the subsystems in the form of the two additional bosonic generators. It would be interesting 
to investigate if the described exotic supersymmetric structure can be related somehow to the 
supersymmetric extensions of the KdV and mKdV equations and corresponding hierarchies that 
are considered within the superspace (superfield) generalizations of the indicated classical (1 + 1)- 
dimensional integrable systems 
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